Product state control of bi-alkali chemical reactions 
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We consider ultracold, chemically reactive scattering collisions of the diatomic molecules KRb. 
When two such molecules collide in an ultracold gas, we find that they are energetically forbidden 
from reacting to form the trimer species K2RD or RD2K, hence can only react via the bond-swapping 
reaction 2KRb — > K2 + Rb2. Moreover, the tiny energy released in this reaction can in principle 
be set to zero by applying electric or microwave fields, implying a means of controlling the available 
reaction channels in a chemical reaction. 

PACS numbers: 34.20.Gj,34.50.Ez,34.50.Lf 



The business of chemical physics is to understand 
the transformation of reactant molecules into product 
molecules during a chemical reaction. On the theory 
side, this is a daunting task, requiring the construction 
of elaborate multidimensional potential energy surfaces, 
complemented by classical, semiclassical, or even fully 
quantum mechanical scattering calculations on these sur- 
faces. For experiments, the goal is to provide as complete 
a picture as possible via the complete determination of all 
initial and final states. For initial states, the use of molec- 
ular beam techniques provides excellent selection of inter- 
nal degrees of freedom in the reactant molecules, and con- 
trol over their relative state of motion. For final states, 
spectroscopic methods can select the relative abundance 
of the different rotational and vibrational states of the 
products, providing a wealth of information from which 
reaction dynamics can be inferred [l|. 

Recently, the experimental attainment of ultracold 
molecules has pushed molecular beam technology to its 
ultimate limit. It is now possible to prepare a molecule 
in a single quantum state in all degrees of freedom, down 
to the nuclear spin state 0, Q . The molecules are more- 
over characterized by an extremely narrow range of ve- 
locities, as set by gaseous temperatures on the order of 
fiK. This circumstance has enabled ordcrs-of-magnitude 
control over chemical reaction rates, by simply altering 
the nuclear spin state of the reactant molecules [1] or else 
by applying an electric field [5[. To probe reactions fur- 
ther, one could imagine transferring the molecules from 
the prepared ground state to any ro-vibrational excited 
state of interest. These advances suggest what is possible 
when one prepares and manipulates the initial states of 
the chemical reaction. 

In this Letter, we suggest that manipulation of the fi- 
nal states may also be possible in certain circumstances. 
Specifically, exit channels can become either energetically 
allowed, or else energetically disallowed, as a function of 
electric field £ . This idea was advanced previously, in the 
context of the H+LiF — > Li + HF reaction, where it was 
argued that the exothermicity of the reaction could be 
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shifted via electric fields 6]. Here we point out that, for 
alkali dimers at ultralow temperature, it is conceivable 
that a previously exothermic reaction can be completely 
turned off. This possibility is afforded by the fact that, 
for alkali dimers, the final states of reaction are not ter- 
ribly different in energy from the reactants. Specifically, 
we consider collisions of a pair of KRb molecules, pre- 
pared in particular rotational states n\ and n 2 . These 
molecules are in general subject to two kinds of reactions: 
the formation of trimers via 

KRb(ni) + KRb(n 2 ) -> K 2 Rb + Rb + A£ trimcr (la) 
KRb(nx) + KRb(n 2 ) -> Rb 2 K + K + A£ trimcr (lb) 

and the bond-swapping reaction 

KRb(m) + KRb(n 2 ) K 2 «) + Rb 2 (n 2 ) + AE hs . (2) 

We report here two circumstances: 1) The reactions 
((T|) that form trimers are energetically disallowed at ul- 
tracold temperatures, as the energy released, A_E t rimor, is 
negative and large - several thousand wave numbers, nine 
orders of magnitude larger than translational kinetic en- 
ergies in the gas. 2) By contrast, the bond-rearrangement 
reactions in produce very small energy differences, 
since the bonds are all covalent and very similar. In zero 
field the reaction ([2]) with n\ — n 2 = is exothermic 
by A-Edimer = 10.4 cm -1 Moreover, the reactants KRb 
are polar whereas the products K 2 and Rb 2 are not, thus 
the products can only be polarized at comparatively high 
fields. Therefore, the relative energy Ai?dimcr(£) between 
reactants and products is, in principle, a function of the 
applied electric field E. Indeed, at fields on the order of 
several xlO 5 V/cm, Ai?b s vanishes and the reaction can 
be turned off altogether. At fields smaller than this, high- 
lying rotational final states can be disallowed, thereby 
changing the possible distribution of product states. 

We begin with the first point, that trimer formation ((T|) 
is energetically forbidden. As the trimer binding energies 
have not been measured, we must calculate them from ah 
initio methods. In general, three spin-1/2 alkali atoms 
can combine to form a trimer with total spin S = 1/2 
(doublet state) or an excited state with S = 3/2 (quartet 
state). Several calculations of doublet states have been 
achieved, for the homonuclear trimers K3 Q and Li3 0] , 
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TABLE I: Molecular properties of the diatomic molecules K 2 , 
KRb, and Rb 2 . Equilibrium bond lengths R e are in A and 
dissociation energies D e axe in cm -1 . (1) denotes a calculation 
performed at the MCSCF+MRCI level while (2) denotes a 
geometry optimization at the RHF-CCSD(T) level accounting 
for BSSE. 



Molecule 


Re (1) 


D e (1) 


Re (2) 


D e (2) 


R e Expt. 


D e Expt. 


K 2 


4.16 


-4293 


3.92 


-4328 


3.92 


-4450.78(15) 


KRb 


4.33 


-4039 


4.05 


-4062 


4.05 


-4217.30(15) 


Rb 2 


4.50 


-3729 


4.18 


-3741 


4.17 


-3993.53(6) 



and for certain molecular Li 2 X systems, where X is an 
alkali- metal atom [9J . Soldan and co-workers have exam- 
ined the two doublet surfaces of K3 at C 2v geometries as 
well as the conical intersection at the equilateral triangle 
geometry described by the D 3 h group Q- 

We employ similar computational techniques as in 
Refs. [7, 9] to compute doublet ground states, within the 
molpro suite of molecular structure codes We use 
the effective core potentials and basis sets of the Stuttgart 
group for the K (ECP10MDF) and Rb (ECP28MDF) 
atoms [ilj . In addition, to adequately model three-body 
forces, we augment these basis functions with a diffuse 
function for each of the s, p, and d orbitals in an even 
tempered manner, as well as a (/-function for the K atom. 

We compare two computational approaches, testing 
them on diatomic alkali molecules. For both approaches, 
we first perform a spin-restricted Hartree-Fock (RHF) 
calculation on the singlet configuration. The first ap- 
proach uses the RHF wave function as the foundation for 
a multi-configuration self consistent field (MCSCF) cal- 
culation [l2l Tl3l| with an active space that includes the 
first excited p orbital of each atom. In addition, we in- 
clude all states resulting from the ns° np 1 configurations 
of each atom. All remaining orbitals are closed, meaning 
they are energy optimized with the restriction that they 
remain doubly occupied. After each MCSCF calculation, 
we perform an internally contracted multi-reference con- 
figuration interaction (MRCI) [3, EH calculation with 
the same active space as in the MCSCF. The minimum 
energy is then compared to the dissociation energy eval- 
uated in the separated-atom limit. 

In the second approach, we perform a coupled 
clusters with single, double and non-iterative triples 
(CCSD(T)) [H[ on the X 1 E state at the minimum R c 
value obtained from the MCSCF+MRCI calculation. We 
then performed a geometry optimization followed by a 
basis set superposition error correction (BSSE)[17| to ex- 
tract the bond length and dissociation energy of the sys- 
tems at the RHF-CCSD(T)+BSSE level of theory. 

To check the adequacy of these computational meth- 
ods, we compare their results with the known binding en- 
ergies of alkali-metal dimers [l8l - [22l |. as shown in Table |U 
These results show that approach 1 gets close to the dis- 
sociation energy of each diatomic species, but noticeably 
overestimates the bond length, especially for the heavier 
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FIG. 1: (Color Online) Contour plots for the doublet ground 
state PES's of K 2 Rb, versus interatomic spacings in an isosce- 
les triangle geometry. The top panel shows the 2 B 2 surface 
while the bottom panel shows the 2 Ai surface. Contours are 
labeled in increments of 750 cm -1 from -5750 - -500 with an 
additional contour near the minimum of each surface. Arrows 
denote the singlet (red) and triplet (green) equilibrium bond 
lengths of the diatomic species. 



alkali-metal systems. By contrast, approach 2 does just 
as well at predicting dissociation energies, but provides 
far better bond lengths. In either method, binding en- 
ergies are clearly reproduced to within several hundred 
cm -1 . We take this as an empirical measure of the cal- 
culations' accuracy. 

We next apply the same methods to determine the 
minima of the three-body potential energy surfaces. The 
minimum of the potential energy surface should lie in the 
(7 2v point group i.e., the odd-atom-out is expected to lie 
on the perpendicular bisector of the line joining the two 
like atoms. Thus we initially restrict ourselves to this 
geometry. As a case study, we will look in detail at the 
K 2 Rb surfaces, but the Rb 2 K surfaces are similar. 

As a first approach we used the RHF-MCSCF+MRCI 
method due to its time efficiency To be consistent with 
the method employed for dimers, we keep the same active 
space, which included the first excited p-orbital of each 
atom. We included 4 states each of 2 B 2 and 2 Ai symme- 
tries, 2 2 Bi, 1 2 A 2 , and 1 4 B 2 states in the MCSCF. In 
each MRCI calculation of the 2 B 2 and 2 Ai surfaces, we 
included 4 states as well as the reference symmetries of 
the other doublets. 

The results of the RHF-MCSCF+MRCI calculations 
are presented in Fig. 1 as contour plots of the poten- 
tial energy surface (PES) in the two independent bond 
lengths i?RbK and Rk 2 , f° r the isosceles triangle geom- 
etry. These surfaces show that the 2 i? 2 surface (upper 
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panel) possesses the global minimum, hence represents 
the ground state of the trimer. The minimum energy oc- 
curs near the singlet bond length of the KRb dimer (de- 
noted by a horizontal red arrow), but is intermediate be- 
tween the singlet and triplet bond lengths of the K2 dimer 
(red and green vertical arrows, respectively). In C2V sym- 
metries, B2 corresponds to an odd reflection about the 
bisector of the isosceles triangle. This means that along 
the K2 bond the electronic wave function should roughly 
resemble the K 2 triplet wave function, and so should the 
the bond length. There is then diminished electron den- 
sity at the center between the two K atoms, and this 
allows the Rb atom to fill this space, bringing the two K 
atoms slightly closer together than they would be in the 
K2 triplet state alone. 

Similarly, the minimum of the 2 Ai surface in the lower 
panel of Fig. 1 is located near the K2 singlet bond length 
(vertical red arrow) , but in between the singlet and triplet 
bonds of the KRb system (horizontal red and green ar- 
rows). The 2 Ai surface requires an even reflection in the 
electronic wave function across the bisector of the isosce- 
les triangle, therefore preferring a singlet-like bond in 
i?K 2 ■ Therefore, the Rb atom does not quite know which 
spin to take since it can form a triplet or singlet with 
one or the other K atoms, but not both. This frustra- 
tion manifests itself with a bond somewhere intermediate 
between singlet and triplet bonds in this coordinate. 

Knowing that the the RHF-CCSD(T)+BSSE results 
in the dimer case are markedly better when compared to 
experimental values, we choose to characterize the min- 
ima of the PES's at RHF-UCCSD(T)+BSSE level of the- 
ory, including a basis set superposition correction. The 
U in UCCSD(T) refers to the spin-unrestricted formal- 
ism due to the open-shell nature of the trimer systems. 
Starting near the parameters obtained above for the 2 B2 
surface, we find the optimum geometry by the method 
of steepest descents. Binding energies are determined at 
the RHF-UCCSD(T) level of theory, accounting for basis 
set overlap errors. The results are presented in Table |TT] 
for all the alkali-metal trimers containing K and Rb. In 
the case of K3, where results have previously been com- 
puted, we find excellent agreement with the reported val- 
ues of the isosceles bond length rj SO = 4.10 A and apex 
angle # apex = 77.13° [7]. Based on these results, we con- 
clude that the formation of the K2Rb trimer would re- 
quire A^trimor = 2(4217) - 5982 = 2452 cm" 1 of energy, 
certainly rendering this reaction impossible at ultralow 
temperatures. 

While a minimum of the surface appears to be in the 
B2 symmetry of the C2 V point group, an examination of 
the trimers K2Rb and Rb2K away from this symmetry 
is needed in order to determine whether this is a global 
minumum. We looked at the systems in the C s point 
group and found minimums which correspond to bent 
geometries using similar optimaization procedures. In 
both trimers the bent geomtry minimum is found to be 
the global minimum making the B2 minimum a saddle 
point. The parameters are given in Table|lTl r^otero refers 



TABLE II: Molecular properties of the alkali-metal triatomic 
molecules containing K and Rb. Bond lengths are in A, an- 
gles are in degrees, and dissociation energies, relative to the 
complete breakup into three free atoms, are in cm -1 . 



Molecule 


Symmetry 


^hctero 


'homo 


e 


D e 


K 3 


C 2 v B 2 


4.10 


5.11 


90.0 


-6059 


K 2 Rb 


C 2v 2 B 2 


4.26 


4.99 


90.0 


-5817 


Rb 2 K 


C 2v ^B 2 


4.26 


5.63 


90.0 


-5533 


K 2 Rb 


C s 


4.25 


4.06 


73.23 


-5982 


Rb 2 K 


C s "A' 


4.24 


4.41 


75.00 


-5611 


Rb 3 


C 2v ^B 2 


4.40 


5.59 


90.0 


-5273 
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FIG. 2: (Color Online) Panel (a) shows the Stark shift for 
the ground and first few rotationally excited states for a pair 
of KRb molecules (colored lines). The black line shows the 
energy of the final states K 2 (n' x = 0) + Rb 2 (n 2 = 0). Panel 
(b) shows the rotational states energetically allowed at zero 
electric field. 



to the shorter of the two bonds between heteronuclear 
pairs when in C s symmetry. 6 is the angle made by A- 
q-B, where A and B are differing atoms and q is the 
midpoint between the homonuclear bond. 

We therefore now consider the bond-swapping reaction 
(|2|) , and in particular the variation of AE\> S with electric 
field. To do so, we look at the energy eigenvalues of the 
reactant and product molecules, as given by the Hamil- 
tonian 



H = B P N 2 - D N A - d ■ £ 



(3) 



where B e is the molecule's rotational constant, D is cen- 
trifugal distortion constant, d is the body-fixed molecular 
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dipole moment (which is zero for the products), and £ 
is the applied electric field. In addition, both kinds of 
molecules also experience a shift due to their electronic 
polarizability a. Using the known KRb dipole moment 
2] and computed molecular polarizabilities (23[, we can 
compute the relative energies of the reactants and prod- 
ucts. 

In Fig. 2a) we plot the energy of the ground state (red) 
and several rotationally excited states of KRb+KRb as a 
function of electric field. Also shown (black) is the ground 
state energy of the products, K 2 (ra = 0) + Rb 2 (n = 0). 
In zero field, the reactants are their natural 10.4 cm -1 
above the products. As the field increases, the energy of 
the polar reactant states decreases rapidly, but that of 
the non-polar product states decreases far more slowly. 
Therefore, at a field beyond ~ 625 kV/cm, the reactants 
are actually lower in energy than the products, Ai?b s 
becomes negative, and the reaction is completely shut 
off. While a static electric field this large is probably im- 
practical to implement, it may be possible to achieve the 
required energy shifts in a suitably designed microwave 
cavity, such as those proposed for trapping polar species 
[2~4| . Theory would then naturally have to account for 
collisions of the field-dressed states [H, HfiJ • 

As a point of reference, Fig. 2b) shows the 235 ener- 
getically allowed states of the products in zero electric 
field, indexed by their rotational quantum numbers n[ 
and n' 2 - It is clear that, before the electric field shuts 



off all reactions completely, it shuts off first higher-n', 
then successively lower-n' states. Recall that the re- 
distribution of the products is one of the key observ- 
ables of physical chemistry. The ability to allow only 
certain values of n' into this distribution will likely pro- 
vide an even greater wealth of information from such ex- 
periments. The detailed effect on chemistry of shutting 
off successively lower rotational exit channels remains to 
be explored. It may be hoped, for example, that the new 
information gleaned would shed additional light on the 
role in these reactions of conical intersections. 

So far we have focused on KRb, since it is the molecule 
for which ultracold chemistry has recently been demon- 
strated experimentally. However, there are other reason- 
able candidates for these experiments, notably RbCs [27j . 
whose reaction is endothermic by 28.7 cm -1 . Such a re- 
action would proceed only by placing the reactant RbCs 
molecules into excited states, which could certainly be 
done. An applied electric field could then still dictate 
which final channels are available, by moving these rota- 
tionally excited states relative to the Rb 2 + Cs2 products. 

Note added in preparation: New work shows that the 
timer reactions (JTJ) are energetically forbidden for all al- 
kali dimer collisions [28j |. 

The authors acknowledge useful discussions with G. 
Quemener and C. H. Greene, as well as funding from the 
NSF. 
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